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^ ■ ABSTRACT 

\ We use the relations between aperture stellar velocity dispersion ((Tap), stellar mass 

■ (Msps), and galaxy size (i?e) for a sample of ~ 10^ early-type galaxies from SDSS/DR7 

to place constraints on the stellar initial mass function (IMF) and dark halo response to 
galaxy formation. We build ACDM based mass models that reproduce, by construc- 
tion, the relations between size, light concentration and stellar mass. Reproducing 
the median dap vs Afsps relation is not possible in models that have both a univer- 
sal IMF and a universal dark halo response. Significant departures from a univer- 
sal IMF and/or dark halo response are required. We show that this degeneracy can 
be broken using the strength of the correlation between residuals of the velocity- 
mass (A log flap) and size-mass (Alogi?e) relations. The slope of this correlation, 
9vR = AlogcTap/AlogiJe, varics systematically with galaxy mass from c)vr — —0.45 
at Msps ~ lO^^M©, to ^vr ^ -0.15 at Msps ^ 10"-^Mq. The virial fundamental 
plane (FP) has cVr = ^1/2, and thus we find the tilt of the FP is mass dependent. 
Reproducing this tilt requires both a non-universal IMF and a non-universal halo re- 
sponse. Our best model has mass-foUows-light at low masses (Afsps ^ IO^-^'^Mq) and 
unmodified NFW haloes at Msps ~ lO^^'^M©. The stellar masses imply a mass depen- 
dent IMF which is "lighter" than Salpeter at low masses and "heavier" than Salpeter 
at high masses. 

Key words: dark matter, galaxies: elliptical and lenticular, cD - galaxies: funda- 
mental parameters - galaxies: haloes - galaxies: kinematics and dynamics - stars: 
luminosity function, mass function 
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1 INTRODUCTION 



The form of the stellar initial mass function (IMF) and the 
response of dark matter haloes to galaxy formation are two 
fundamental unknowns, which are important in many areas 
of astrophysics. For example, the IMF is needed in order to 
convert observations of integrated stellar light into stellar 
masses and star formation rates (two fundamental parame- 
ters in galaxy evolution studies), and for modelling the su- 
pernova rates, chemical evolution, and production of ioniz- 
ing photons in galaxies. The response of dark matter haloes 
to galaxy formation is needed in order to constrain the na- 
ture of dark matter from observations of the structure of 
dark matter haloes (e.g., central density slopes, central den- 
sities) and indirect dark matter detection experiments. Halo 
response also plays a key role in determining the fundamen- 
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tal scaling relations between velocity and light of galaxies 
(Faber & Jackson 1976, TuUy & Fisher 1977) 

For many years the IMF and dark halo response were 
thought to be universal: The IMFs of external galaxies 
are the same as measured in the Milky Way (Kroupa 
2001, Chabrier 2003); and dark matter haloes contract adi- 
abatically in response to galaxy formation (Blumenthal 
et al. 1986). However, recent observations and numerical 
simulations have cast doubt on these assumptions. We note 
that while there have been claims of observational evi- 
dence for dark halo contraction in early-type galaxies (e.g., 
Schulz et al. 2010; Napolitano et al. 2010; Trujillo-Gomez 
et al. 2011; Chae et al. 2012) these conclusions depend (triv- 
ially) on the assumption of a universal Milky Way type IMF 
(e.g.. Button et al. 2011a). 

It has been known for many years that baryonic ef- 
fects can, in principle, result in reduced halo contraction or 
even halo expansion (e.g., Navarro, Eke, & Frenk 1996; El- 



© 2012 RAS 



2 Button et al. 



Zant et al. 2001). But only recently have these effects been 
demonstrated in fully cosmological simulations of galaxy for- 
mation (Johansson et al. 2009; Abadi et al. 2010; Duffy 
et al. 2010; Governato et al. 2010; Maccio et al. 2012; Mar- 
tizzi et al. 2012). Although some authors maintain that dark 
halo contraction is an unavoidable consequence of galaxy 
formation in ACDM cosmologies (e.g., Gnedin et al. 2011). 

The traditional route to constraining the IMF is 
through direct star counts. In the Milky Way this has 
shown the IMF to have a power-law shape dN/dm oc m~^, 
with X ~ —2.3 at masses above m ~ IMq (Salpeter 
1955), and a turn over at lower masses (Kroupa 2001; 
Chabrier 2003). Since most of the mass is in low-mass stars, 
Kroupa/Chabrier IMFs have lower stellar mass-to-light ra- 
tios (by about 0.2-0.25 dex) than a Salpeter IMF. Outside 
of our Galaxy counting individual stars is usually not fea- 
sible (especially down to the low masses required to fully 
constrain the IMF) . There are several approaches which can 
be used to probe the IMF at different mass scales in extra- 
galactic systems. 

Constraints on the shape of the IMF at high masses 
typically involve measurements of Ha and/or far ultraviolet 
(FUV) fluxes. For example, the equivalent width of Ha can 
constrain the slope of the IMF above ~ IMq (Kennicutt 
1983, Hoversten & Glazebrook 2008), and the ratio between 
Hq and FUV fluxes probes IMF above ~ lOM© (Meurer 
et al. 2009; Lee et al. 2009). Broadly speaking these meth- 
ods show that the high mass IMF in external galaxies is 
similar to that of the Milky Way, but there are systematic 
discrepancies which could be explained by variation in the 
slope and/or upper mass limit of the IMF. Note that while 
modest changes in the upper end of the IMF don't impact 
the stellar mass-to-light ratios of old stellar populations, top 
heavy IMFs can leave behind enough stellar remnants to sig- 
nificantly increase stellar mass-to-light ratios. 

The low mass end (below ~ IMq) of the IMF can be 
probed by comparing the mass (or an upper limit) of all the 
stars from dynamics and/or strong gravitational lensing to 
the stellar mass obtained from stellar population synthesis 
(SPS) models. Such studies require IMFs lighter (i.e., lower 
stellar mass-to-light ratios) than Salpeter in spiral galax- 
ies (BeU & de Jong 2001; Bershady et al. 2011; Button 
et al. 2011a, Barnabe et al. 2012; Brewer et al. 2012), but 
are consistent with a universal Milky Way type IMF. 

In elliptical and lenticular galaxies, Cappellari 
et al. (2006) showed that a Salpeter IMF over-predicts 
the dynamical mass-to-light ratios of some "fast-rotating" 
galaxies. Using dynamics and strong gravitational lens- 
ing from the SLAGS survey (Bolton et al. 2006) Treu 
et al. (2010) showed that in fact a Salpeter IMF is allowed 
in massive (stellar mass Matar ^ 10^^ Mq, velocity disper- 
sion a > 200kms~^) elliptical galaxies, even accounting 
for "standard" dark matter haloes. Extending this study 
to include constraints from weak lensing, predictions from 
AGDM , and the possibility of dark halo contraction. Auger 
et al. (2010b) concluded that a Salpeter-type IMF was 
strongly favored over Milky Way type IMF. In order to 
minimize uncertainties in subtracting off the dark matter. 
Button, Mendel, & Simard (2012) studied the densest 
~ 3% of early-type galaxies in the SDSS. These galaxies 
have a fundamental plane correlation consistent with no 
dark matter within an effective radius. The average IMF 



of these galaxies is close to Salpeter, with evidence that 
redder/bluer galaxies have heavier/lighter IMFs. 

Another method for probing the low mass end of the 
IMF (in non star forming galaxies) is by measuring the 
strength of dwarf star sensitive absorption lines. Apply- 
ing this method to nearby massive elliptical galaxies van 
Bokkum & Gonroy (2010) found an IMF steeper (below 
~ IMq) than that of the Milky Way, and most likely steeper 
than Salpeter. 

Thus there is good agreement from different methods 
that the IMFs of the most massive or most dense elliptical 
galaxies are heavier than that of the Milky Way, and most 
likely similar to a Salpeter IMF below ~ IMq . The next step 
is to constrain how the IMF varies in a more general sam- 
ple of early-type galaxies. The ATLAS3B project has made 
progress towards this aim (Gappellari et al. 2012), but their 
study is limited by the inherent degeneracies between the 
baryonic and dark matter components in the mass model- 
ing of galaxies (e.g.. Button et al. 2005). Specifically there 
is a factor of ~ 2 range in derived stellar mass-to-light ra- 
tios between models in which mass-foUows-light and models 
with cosmologically motivated adiabatically contracted dark 
matter haloes (see also Button et al. 2011a). 

In this paper we place constraints on the mass depen- 
dence of the IMF and dark halo response using scaling rela- 
tions of a large sample (> 10^) of early-type galaxies from 
the Sloan Bigital Sky Survey (SBSS, York et al. 2000). 
This approach has strengths and weaknesses compared to 
detailed mass models of individual galaxies (e.g.. Button 
et al. 2011b; Sonnenfeld et al. 2012; Gappellari et al. 2012). 
The weakness is we only get constraints for average galax- 
ies (of a given mass), the strength is that our constraints 
are free of statistical uncertainties (in e.g., stellar mass-to- 
light ratios, dark halo masses, dark halo concentrations, etc) 
that necessarily effect studies with smaller samples. Both ap- 
proaches are valid, providing complementary constraints to 
the variability of the IMF and dark halo response. 

We construct dynamical models that consist of spheri- 
cal distributions of stars and dark matter. These models are 
constrained to reproduce a number of observational and the- 
oretical constraints (e.g.. Button et al. 2011a). These models 
have two unknowns: the stellar IMF and the response of the 
dark matter halo to galaxy formation. We use the observed 
stellar velocity dispersion - stellar mass (Faber-Jackson) re- 
lation to single out allowable combinations of IMF and halo 
response. We show that the IMF and/or the halo response 
has to vary with galaxy mass. However, from the velocity 
dispersion vs stellar mass relation alone, we are unable to 
uniquely determine which if either should vary, or the abso- 
lute normalization of the IMF and halo response. 

Our additional constraint is the slope of the correlation 
between the residuals of the velocity - stellar mass (A log a) 
and size - stellar mass (Alogi?) relations, which for brevity 
we denote 9vr = A log a/ A log 7?. This is a useful con- 
straint because it depends on the dark matter fraction within 
(roughly) the half-light radius (Gourteau & Rix 1999; But- 
ton et al. 2007). If there is no dark matter then Ova = —1/2, 
while dark matter dominated galaxies are expected to have 
dvK > 0. As we show below dvn is related to a more famil- 
iar concept - namely the tilt of the fundamental plane (FP, 
Bressler et al. 1987; Bjorgovski & Bavis 1987). 

This paper is organized as follows. The relation between 
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dvR and the tilt of the FP is shown in §2, the main observed 
scahng relations are presented in §3, with additional details 
given in appendix A. Constraints on the IMF and halo re- 
sponse from the velocity - mass relation and the tilt of the 
FP are discussed in §4, and §5 respectively. A summary is 
given in §6. Appendix A gives the details of our mass mod- 
els. We adopt a ACDM cosmology with f^A = 0.7, JIm = 0.3 
and Ho = 70 kms"^ Mpc"^ 



2 THE TILT OF THE FUNDAMENTAL PLANE 

The fundamental plane is a well-known scaling relation be- 
tween the size (R), and surface brightness (/), and velocity 
dispersion (a) of early-type galaxies: 

Rxa"!''. (1) 

For r-band sizes and surface brightnesses, the exponents are 
(a, 6) = (1.49±0.05, -0.75±0.01) (e.g. Bernardi et al. 2003). 
This differs from the simplest mass-foUows-light model, 
known as the virial plane, in which (a, &) = (2,-1). The 
differences in exponent are referred to as the tilt of the fun- 
damental plane. 

The origin of the tilt has been the subject of debate 
for many years (e.g., van der Marel 1991; Bender, Burstein 
& Faber 1992; Ciotti, Lanzoni, & Renzini 1996; Graham & 
CoUess 1997; Prugniel & Simien 1996, 1997; Pahre, Djor- 
govski, & de Carvalho 1998; Padmanabhan et al. 2004; Tru- 
jillo, Burkert, & BeU 2004; Robertson et al. 2006; Bohon 
et al. 2007, 2008b; Graves & Faber 2010). There are two 
principle explanations for the tilt: (1) Variation of total mass 
to light ratio; and (2) Structural non-homology. Using the 
fact that / oc L/R^ one can re- write Eq. [T]as 

Ro^/Lo^o'^T^-''. (2) 

Thus the virial plane corresponds to Ra^ /L — constant, and 
the tilt of the fundamental plane corresponds to variation in 
Ra'^/L. 

In the case of structural homology the total mass, Mtot, 
within radius R, is directly proportional to Rcr^ , and thus 
variation in Ra^ /L corresponds to variation in Mtot/ L. The 
total mass-to-light ratio can vary as a result of variation in 
stellar mass-to-light {Msta,i/L) or variation in the total mass 
to stellar mass (Mtot/Afatar) ratio. The former is a result of 
stellar population differences, while the latter is a result of 
variation of dark matter fractions. 

In the case of constant Mtot/L, structural non- 
homology results in the total mass not being directly pro- 
portional to Ra^ . Examples of structural non-homology are 
variation of the stellar mass profile (e.g., a variation in the 
Sersic index or concentration parameter), and variation in 
the anisotropy of stellar orbits. 

Going back to the fundamental plane as written Eq. [T] 
both coefficients are tilted with respect to the virial predic- 
tion. It turns out that the fundamental plane can be written 
so that the tilt is in a single coefficient. We can re- write Eq.[T] 
as 

o oc L'ii'* (3) 

where c = ~b/a, and d = (1 -|- 2b) /a. The correspond- 
ing virial coefficients are (c, d) = (1/2,-1/2). The expo- 
nents from Bernardi et al. (2003) result in (c, d) = (0.50 ± 



0.02, — 0.33±0.01), and thus the tilt of the r-band fundamen- 
tal plane is consistent with being entirely due to d 7^ ^1/2- 
Replacing luminosity with stellar mass in Eq. O we see that 
the strength of the correlation between the velocity-mass 
and size-mass relations is equivalent to the tilt of the fun- 
damental plane, i.e., 9vr = d. 



3 OBSERVED SCALING RELATIONS 
3.1 Sample Overview 

Here we give an overview of our observational sample of 
early-type galaxies. We use a similar (but not identical) 
selection procedure as described in more detail in Button 
et al. (2011a). In summary we use four selection criteria: (1) 
A spectroscopic redshift 0.005 < z < 0.3; (2) A spectrum 
classified as early-type by SDSS (eCLASS < 0); (3) A red 
{g — r) color, based on valley in color - stellar mass plane; 
(4) A minor-to-major axis ratio, h/a > 0.5 (which removes 
dusty edge-on spirals) . These cuts are designed to select non- 
star forming galaxies, so our sample includes ellipticals and 
lenticulars. We also apply a redshift dependent minimum 
stellar mass to remove the color bias at fixed stellar mass 
which results from the SDSS r-band magnitude limit for 
spectroscopy. Our final sample consists of ~ 150 000 early- 
type galaxies. 

There are two velocity dispersion measurements pub- 
licly available for SDSS galaxies. We refer to these as "SDSS" 
and "Princeton". At most masses these two measurements 
are in excellent agreement, in addition the scatter in these 
two measurements is in good agreement with the reported 
uncertainties. However, at low masses (and low dispersions) 
there are significant differences (see also Hyde & Bernardi 
2009). In what follows, for the observed velocity dispersion 
we adopt the logarithmic average of the two measurements. 
Where appropriate we will use the individual measurements 
to gauge systematic uncertainties in scaling relations. 

When constructing the Faber- Jackson relation and the 
fundamental plane from fiber spectroscopy it is customary 
to correct the aperture velocity dispersions to some fiducial 
radius, such as the effective radius {Re), or one eighth of 
the effective radius. However, doing so requires making an 
assumption for the variation of o with radius. Since we will 
directly use the strength of the correlation between a and 
Rc at fixed mass to constrain our models, we do not wish to 
introduce an artificial correlation. Instead, our approach is 
to measure observed scaling relations using the uncorrected 
fiber velocity dispersions, (Tap, and to explicitly compute the 
velocity dispersion within the SDSS fiber in our models using 
the spherical Jeans equations. 

We use stellar masses from the MPA/JHU databasfEl- 
These are derived by fitting ugriz photometry with Bruzual 
& Chariot (2003) SPS models assuming a Chabrier (2003) 
IMF. To account for a non-universal IMF, or systematic er- 
rors in stellar mass measurements, we allow for an offset 
between the true stellar masses, Mstar and these SPS stellar 
masses, Msps, which we denote Aimf = log(Mstar/Msps). 

The uncertainty on the stellar masses are nominally 
~ 0.1 dex, but there will be some galaxies with much larger 

^ Available at |http://www.mpa-garchlng.mpg.de/SDSS/DR7/| 
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Figure 1. Observed velocity - stellar mass (VM, left) and size - stellar mass (RM, right) relations for SDSS early-type galaxies. Velocity 
dispersions (cap) are measured within the SDSS 3 arcsec diameter aperture, the stellar masses (Mgps) assume a Chabrier (2003) IMF, 
the sizes (-Re) are circularized half-light radii measured in the i — band. For the VM relation the black solid lines show the relations 
obtained for the two different velocity dispersion measurements, which are only visibly different at low masses. Upper panels: The solid 
points show medians in bins of 0.2 dex width in Mgps. The solid red lines show double power-law fits to these points (parameters in 
Table IX2t . while the dotted lines show the 16th and 84th percentiles. Lower panel: The open points show the scatter in each mass bin. 
The red dotted lines show fits to these points (parameters given in Table IX2t . 



uncertainties. At the highest masses, and largest size-ofFsets, 
these outliers can bieis the slopes of the scaling relations. 
As in Button, Mendel, & Simard (2012) we clean the sam- 
ple by using the relation between stellar mass-to-light ratio 
and velocity dispersion: Msps/ir vs (Tap. We remove ~ 700 
galaxies (0.5% of our ETG sample) that are more than ±4cr 
from the median relation. This number is ~ 100 times higher 
than the expected number of high-a offsets for log-normal 
scatter in Msps/Ir- 



3.2 Mean Relations 

The observed velocity dispersion - stellar mass (VM) and 
size - stellar mass (RM) relations for our sample of SDSS 
early-type galaxies are shown in Fig. [1] The upper panel 
shows the median (filled circles) together with the 16th and 
84th percentiles (dotted lines) of velocity and size in bins 
of stellar mass. The median relations are well fitted with 
double power- laws (in the variables 10^ and 10^): 



X /3 - a , 
y = yo + oi{x - a::o) H log 

7 



4- -!-lo^'^~^i^ 



(4) 



Where a is the slope at a; ^ xo, /3 is the slope at a; ^ a;o, 
xo is the transition scale, j/o = ^(2:0), and 7 controls the 
sharpness of the transition (higher 7 results in a sharper 
transition). The parameters of the fits are given in Table IX2I 
For the velocity-mass relation the slope varies from ~ 0.4 at 
low masses to ~ 0.2 at high masses, while for the size-mass 
relation the slope varies from ~ at low masses to ~ 0.7 at 
high masses. 



The lower panels show the scatter about the median re- 
lations (open circles). For both velocity-mass and size-mass 
relations the scatter decreases from low to high mass. The 
mass dependence of the scatter is fitted with the following 
function: 



y 



2/2 + ■ 



(5) 



1 + 10^(^ 

Here j/i is the asymptotic value at x ^ a;o, J/2 is the asymp- 
totic value at X ^ xo, xo is the transition scale, and 7 con- 
trols the sharpness of the transition. The values of the best 
fit parameters are given in Table IX2I 

3.3 Residual Relations 

The correlation between the residuals of the VM and RM 
relations, A log i?c and A log aap, for the full sample of early- 
type galaxies is shown in Fig. [2] The histograms show the 
distribution of Alog_Re and AlogtJap. The standard de- 
viation of the velocity and size residuals is 0.08 dex and 
0.13 dex, respectively. The filled circles show the median 
A log (Tap in bins of Alogi?c, while the dotted lines show 
the 16th and 84th percentiles of the distribution. The me- 
dian errors on velocity, size, and mass are indicated with 
the compass, and shows that most of the scatter about the 
median relation is due to observational errors. 

A power-law fit to these data results in a slope of 
i9vR = —0.29 ± 0.01. The virial fundamental plane has 
i9vR = —0.5, which is shown with the diagonal dashed line. 
Thus the fact that we observed 9vr 7^ —0.5 is equivalent 
to the statement that the observed fundamental plane is 



© 2012 RAS, MNRAS OOO.ITlfTil 




\T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 
\ 

" \ ,9.4 - 1 1.8 slope=-0.29 




A log,o [Re/(kpc)] I M 

SPS 

Figure 2. Correlation between residuals of the velocity-mass and 
size-mass relations (Fig. [TJ for all early-type galaxies. The slope 
of the best fit linear relation is cJyR = —0.29. The long-dashed 
line is the virial relation with slope= —0.5. The compass shows 
the direction the median errors on velocity (o"ap)i size (Re), and 
mass (Mgps) scatters galaxies in this plane. This shows that most 
of the scatter about the median relation is due to observational 
errors. 

tilted with respect to the virial fundamental plane. The slope 
varies with Alog J?e, indicating that the tilt of the funda- 
mental plane is not a constant. For large negative size off- 
sets, i9vR — —0.5, suggesting that baryons dominate within 
the effective radius, whereas for large positive size offsets, 
9vR — 0, suggesting that there is significant dark matter 
(Courteau & Rix 1999). 

Fig.[3]shows dvn computed in stellar mass bins of width 
0.2 dex. This shows that while dvR is always negative, it is 
not a constant, i.e., the tilt of the fundamental plane in- 
creases with increasing stellar mass. As an estimate of sys- 
tematic uncertainties we measure dvn using the two sets 
of velocity dispersions available for SDSS galaxies. The two 
measurements are shown with the solid lines in Fig. [S] They 
are in good agreement above Msps ~ lO^^M©, but below 
this mass there are significant differences. 

As a final note we stress that the dyn shown in Fig. [3] 
are not corrected for aperture effects or measurement errors. 
Our approach is to explicitly include these effects in our 
models, rather than try to correct for them in the data. We 
will show that an observed dvR ~ —0.3 can be explained by 
models in which mass follows hght (and thus follow the viral 
fundamental plane between a^, R^, and A/star). 



4 CONSTRAINTS FROM THE 
VELOCITY-MASS RELATION 

In this section we use the velocity-stellar mass relation to 
constrain the two free parameters of our model (see Button 



logio [MsPs/(M0)] 

Figure 3. Slopes of the correlation between residuals of the 
velocity-mass (A log (Tap) and size-mass (Alogi?e) relations, as 
a function of stellar mass. The horizontal bar shows the width 
of the bin in stellar mass, while the vertical error bar shows the 
uncertainty on the slope. The two solid lines show the results ob- 
tained using the two different set of stellar velocity dispersions, 
indicating systematic errors at low masses. The virial FP has a 
slope of -0.5, and thus we find that the tilt of the FP is mass 
dependent. 

et al. 2011a and Appendix A for more details); the stellar 
mass normalization Aimf ~ log(Mstar/Msps) and dark halo 
response v. 

4.1 Universal IMF and universal halo response 

We start by constructing model samples with a universal 
Chabrier IMF. We consider four halo responses: standard 
adiabatic contraction f = 1 (Blumenthal et al. 1986), re- 
duced halo contraction u — 0.5 (c.f., Abadi et al. 2010); 
no halo contraction v = (i.e., NFW haloes); and halo 
expansion v = —0.5. In addition we consider a model in 
which mass-foUows-light (MFL). For each model we com- 
pute aperture velocity dispersions for 5000 model galaxies 
evenly spaced in log(Afsps/AfQ) (from 9.3 to 11.9), and in- 
cluding log-normal scatter in sizes, dark halo masses, and 
dark halo concentrations. We then re-sample these galax- 
ies (100 times) according to the observed intrinsic distribu- 
tion of stellar masses, and add measurement errors in stellar 
mass, size and velocity dispersion (see Appendix A). This 
procedure results in a sample of ~ 150000 model galaxies 
with the same distribution of stellar masses and sizes as our 
observed sample. 

The median velocity-mass relations of these models are 
shown in Fig. |4l None of our models is able to reproduce 
the observed VM relation, even allowing for zero point off- 
sets (corresponding to different, but still universal IMFs). 
There are two primary solutions to this problem: 1) Allow 
the stellar mass to vary from that predicted by a univer- 
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"oQio (Msps/Mq) logio (Msps/M©) 

Figure 5. Left: VM relations for data (blaclc filled circles) and mass-foUows-light (MFL) models (red and blue lines) assuming a Chabrier 
(2003) IMF. The model galaxies have constant velocity anisotropy with (3 = (red solid line) and /3 = 0.5 (blue dashed line). Right: 
Relation between stellar mass (Afstar) required for the model to match the observed VM relation and stellar population synthesis mass 
(Mgps). The offset between these two masses is well fitted with a power-law, with parameters as indicated. Interpreting this offset in 
terms of the IMF requires "heavier" IMFs in more massive galaxies. 
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Figure 4. Comparison between observed and model VM rela- 
tions assuming a Chabrier IMF. The filled circles and solid black 
line shows the observed median relation, with error bars indicat- 
ing the error on the median. The colored lines show the model 
relations assuming isotropic stellar velocity dispersions (/3 = 0). 
The models only differ in the structure of the dark matter halo, 
which ranges from adiabatically contracted NFW dark matter 
haloes {y = 1, blue line) to no dark matter (MFL, magenta line). 
None of these models reproduces the slope of the VM relation, 
indicating the need for a non-universal IMF and / or non-universal 
dark halo response. 



sal Chabrier IMF, with "heavier" IMFs in more massive 
(i.e., higher Msps) galaxies; or 2) Allow the halo response 
to vary with galaxy mass, with e.g., no halo contraction in 
low mass galaxies and stronger halo contraction in progres- 
sively higher mass galaxies. 

In principle, another solution would be to allow the stel- 
lar anisotropy to vary with galaxy mass, however, observa- 
tions find no evidence for a mass dependence to the stellar 
anisotropy (e.g., Gerhard et al. 2001). Furthermore, as we 
show below, the maximum range of reasonable anisotropy 
only has a ~ 10% eflect on the derived masses. 

We now construct models with non-universal IMFs 
and/or non-universal halo responses that reproduce the me- 
dian VM relation. In § Owe will use the tilt of the funda- 
mental plane (Fig. |3]) to distinguish between these models. 



4.2 Non-universal IMF with mass-follows-light 

For MFL models we can easily calculate the stellar masses 
required to match the observed velocity dispersions using 



Afel 



MsPs(o"ap/(Tap,SPs) 



(6) 



where CTap is the observed velocity dispersion and (Tap.SPS 
is the model velocity dispersion computed assuming Msps- 
The stellar mass offset, Aimf, as a function of SPS mass for 
two MFL models (with different anisotropy) is shown in the 
right panels of Fig. [5] For 10^'^ < A'/sps < lO"'^ the stellar 
mass offset is well fitted with a linear relation: 

Aimf = log(A/,tar /A/sps) = a -f 61og(AfsPs/10"A//Q). (7) 

For both isotropic (/? = 0) and radially anisotropic orbits 
(with j3 = 0.5) we find a slope of 6 ~ 0.14. This results in 
a factor of ~ 2 difference in Aimf across the range of Afsps 
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Figure 6. Comparison between observed (black circles) and model (red lines) VM relations. Each model has a different halo response 
(parametrized by u) and stellar mass normalization (parametrized by a and b). Models with stronger halo contraction require lighter 
IMFs (lower a). For a universal halo response (i.e., fixed u) a non- universal IMF (b 7^ 0) is required. 



Table 1. Parameters of models. 



Halo Response /3 Aimf zero Ajmp slope 



V 




a 


b 


Variable 


0.0 


0.08 


0.00 


1.0 


0.0 


0.03 


0.14 


0.5 


0.0 


0.15 


0.14 


0.0 


0.0 


0.22 


0.14 


-0.5 


0.0 


0.26 


0.14 


MFL 


0.0 


0.28 


0.14 



that we study. In terms of normalization, for /? = (0.0, 0.5) 
we find a = (0.278,0.238), which is slightly higher than a 
Salpeter IMF (Aimf 0.23 ± 0.01). 

4.3 Degeneracy between IMF and halo response 

For models with dark matter haloes the model velocity dis- 
persion depends non-trivially on the stellar mass, and thus 
when we change Aimf or the halo response, we need to re- 
calculate the model. We determine the IMF offset param- 
eters (a, 6) for such models with a grid search. The stel- 
lar mass offset parameters for our suite of models are given 
in Table [1] and the VM relations are shown in Fig. |6] All 



models require a variable IMF, but the mass dependence, 
h ~ 0.14, is consistent with being independent of the halo 
response model. At the reference stellar mass (Msps = 
10^^ M0) standard adiabatic contraction {u = 1, Blumen- 
thal et al. 1986) requires an IMF close to Chabrier, while 
no dark halo contraction {v = 0) requires an IMF close to 
Salpeter. This is in agreement with our earlier study Button 
et al. (2011a). 

One can construct models that match the VM relation 
with a universal IMF, but a non-universal halo response. An 
example of such a model is shown in the lower right panel of 
Fig. [S] The IMF normalization a = 0.08 is slightly heavier 
than Chabrier (2003) . The halo response parameter depends 
on stellar mass as 1/ = 1.0 + 1.0(log Msps — 11), with a max- 
imum of f = 1. So this model has adiabatic contraction for 
Msps > lO^M©, and no-contraction at Msps = 10^°Mq. 

The most relevant study in the literature to ours 
is by Auger et al. (2010a). These authors fitted ACDM 
based models to a sample of ~ 50 massive (Msps = 
^gii.35±o.20^^'j giiiptical galaxy strong lenses from the 
SLAGS survey (Bolton et al. 2006). Observational con- 
straints were masses from strong lensing, aperture stellar 
velocity dispersions from SDSS, and weak lensing data from 
HST. Auger et al. (2010a) parametrized the relation between 
Msps and Mstar by logMsps = logM^tar - r7(logMstar - 
11) — a. In terms of our parametrization: a = a/{l — rj) and 
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Figure 7. Slopes of the relation between residuals of the VM and RM relations {c?vR,) vs stellar mass (normalized to a Chabrier IMF) 
for models in Fig. [G] The data (black points) is as in Fig. [S] For the models the shaded region corresponds to the fiducial measurement 
errors, while the dashed region shows 20% higher and lower errors on stellar masses. For all but the highest masses the mass-foUows-light 
(MFL) model (magenta, upper left panel) fits the data the best. For the highest masses, models with no halo contraction are favored, 
but mild halo contraction is consistent with the data. Thus reproducing the tilt of the fundamental plane requires both and non-universal 
IMF and non-universal halo response. 



b = r]/(l — rf). Their best fitting model liad uncontracted 
NFW iialoes and a mass dependent IMF witli a = 0.03±0.03 
and Tj — 0.08 ± 0.04, wiiere tlie reference IMF was Salpeter 
(1955). Converting to a Chabrier (2003) IMF results in 
a = 0.26 ± 0.06, b = 0.09 ± 0.06, whicli is consistent with 
this paper for the masses our respective studies overlap. 



5 CONSTRAINTS FROM THE TILT OF THE 
FUNDAMENTAL PLANE 

We now use the tilt of the fundamental plane to break the 
degeneracy between IMF and halo response. Fig. [7] shows a 
comparison between Ova and stellar mass for models and 
observations. The observations are shown with black points, 
lines and error bars (as in Fig. [31). The models are shown with 
colored shaded regions, as indicated. The models include 
measurement errors in stellar masses, sizes and velocity dis- 
persions. The effect of measurement errors is to increase the 
strength of fundamental plane tilt (i.e., i9vr becomes less 
negative). Errors on stellar masses have the strongest effect 
on dvK- To indicate the dependence of our results to stellar 
mass errors the dashed lines show models with 20% higher 
and lower measurement errors on stellar masses. 



The first thing that is apparent is that all of the models 
predict that Ovr should vary with stellar mass, with (gen- 
erally) more negative dvn at lower masses, in qualitative 
agreement with the observations. However, none of the mod- 
els reproduces the shape and normalization of (?vr vs Msps 
in detail. In particular models with standard adiabatic con- 
traction (blue line, lower middle panel) or a universal Milky 
Way type IMF (grey line, lower right panel) are the most 
discrepant, under predicting the strength of Ovr by up to 
~ 0.2. These models have the "lightest" IMF normalizations, 
and hence the highest dark matter fractions (~ 50%) within 
an effective radius. The increased dark matter has the effect 
of dampening the changes in velocity that are expected from 
changing the sizes of the galaxies. 

Below a stellar mass of Msps ~ lO"'^ the MFL model 
(magenta line, upper left panel) is the only one that is con- 
sistent with the data. At the lowest masses (Msps ~ 10^'^) 
there is significant uncertainty in the observations, and 
while MFL models fit best, they are not required. Above 
Msps ~ 10^^'^ the MFL model progressively over-predicts 
the strength of dvn- Reproducing the observations thus re- 
quires a mass dependent halo response: expansion (red line) 
for masses lO"'^ < Msps < lO"", and uncontracted NFW 
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Figure 8. Summary of the main results for the steUar IMF and halo response of this paper. Reproducing the average Velocity - Mass 
is degenerate between IMF and dark halo response. For a universal dark halo response a non-universal IMF is required (dotted blue 
lines), with Mgtar/AfsPS -'^^SPS ' "heavier" IMFs in higher mass galaxies. The tilt of the stellar mass fundamental plane enables 
us to break this degeneracy. Our best model (thick black line) has mass follows light at low masses Mspg < 10^^'^, expansion for 
Mgpg ~ iQii-^"!! '' and uncontracted NFW haloes for Mgpg ~ 10^^ ®. At the highest masses A/gpg > 10^^ * halo contraction is allowed 
given reasonable uncertainties in stellar mass errors. Changing the stellar anisotropy from /3 = (isotropic orbits) to /3 = 0.5 only 
lowers the derived stellar masses by ~ 10%, as indicated by the error bar in the top left corner. The cyan dot shows the result for dense 
early-type galaxies from Dutton et al. (2012), while the green error bar shows the results for strong lenses from Auger et al. (2010a). 
The upper horizontal axis shows the average aperture stellar velocity dispersions at the corresponding stellar masses. 



haloes (green line) for Mgps ~ lO'^^'^. Given reasonable un- 
certainties in the stellar mass errors, models with contracted 
dark matter haloes (cyan and blue) are consistent with the 
data in the most massive galaxies Mgps ^ lO^^ *. 

The fact that MFL does not reproduce the fundamental 
plane for the most massive galaxies is consistent with results 
from strong gravitational lensing, which find average total 
mass density slopes that are close to isothermal 7 ~ — 2.08± 
0.03 (Auger et al. 2010b), and thus shallower than MFL 
(7 ~ —2.3 for a Hernquist profile). 

It might seem surprising that MFL models, which fol- 
low the virial fundamental plane, can reproduce the observed 
fundamental plane tilt. This can be explained by two effects. 
The first, as mentioned above, is that measurement errors 
increase the observed fundamental plane tilt (i.e., they make 
dvR less negative). The second is a result of measuring the 



velocity dispersions within a fixed physical aperture, rather 
than within a relative aperture such as an effective radius. 
Each of these effects accounts for roughly half of the ob- 
served tilt of the fundamental plane for intermediate mass 
galaxies. 



6 SUMMARY 

We use the relations between aperture stellar velocity disper- 
sion (o"ap), and stellar mass (Msps), and galaxy size {Rc) to 
constrain the dark halo response and stellar IMF in early- 
type galaxies. Our observational sample consists of ~ 150 
000 early-type galaxies from the SDSS/DR7 at redshift z ~ 
0.1. We construct mass models with ACDM haloes (modified 
by contraction or expansion), as well as mass-foUows-light 



© 2012 RAS, MNRAS 000.[TlfT4l 



10 Button et al. 



(MFL) , that reproduce the observed distribution of SPS stel- 
lar masses and galaxy sizes. Our models include uncorrelated 
scatter in galaxy sizes, halo masses and halo concentrations. 
The remaining free parameters of the model are constrained 
by the median relation between (Jap and Msps (the Faber- 
Jackson relation), and the strength of the correlation be- 
tween residuals of the (Tap vs Msps and effective radius i?c 
vs Msps relations: 9vr = A log cr ap/ A log _Ro (equivalent to 
the tilt of the fundamental plane) . We summarize our results 
in Fig. [8] and as follows: 



• Reproducing the median (Tap — Msps relation is not pos- 
sible for models with both a universal IMF and universal dark 
halo response. Significant departures in either/or both are 
required. 

• Models with a universal halo response to galaxy forma- 
tion (as well as MFL) require heavier IMFs in higher mass 
galaxies and are consistent with Mstar oc Msps ■ ^ stellar 
mass of Msps = W^^ Mq these models with adiabatic con- 
traction require close to Chabrier IMFs, models with uncon- 
tracted NFW haloes require close to Salpeter IMFs, while 
MFL models require IMFs heavier than Salpeter. 

• A model with a universal (close to Kroupa) IMF and 
mass dependent dark halo response ranging from expansion 
at low masses (Msps <■ 10^" Mq) to adiabatic contraction 
at high masses (Msps <; 10^^ Mq) is able to reproduce the 
(Tap — Msps relation. 

• We find that dvn varies systematically with stellar 
mass. The minimum is dvR ~ —0.45 at A/sps = W^'^Mq 
(Chabrier IMF) and increases to dvR — —0.15 at Msps = 
IO^^ ^Mq. The virial fundamental plane has dyn = —1/2. 
And thus the tilt of the fundamental plane is not a constant, 
as is generally assumed. 

• The MFL model successfully reproduces the mass 
dependent fundamental plane tilt for masses Msps < 
1O"-^M0. However, above Msps ~ 10"'^ the MFL model 
progressively over-predicts the strength of dvn. Reproduc- 
ing the observations favors expansion for masses 

Msps ^ 10"", and uncontracted NFW haloes for Msps > 
10^^'^. At the highest masses models with halo contraction 
are also consistent with the data. 

• Models with a universal IMF (and non-universal halo 
response) are unable to reproduce the tilt of the fundamental 
plane. 



Our results are in agreement with several recent stud- 
ies which also favor IMFs that are "heavier" than Salpeter 
in massive elliptical galaxies (van Dokkum & Conroy 2010; 
Auger et al. 2010; Sonnenfeld et al. 2012). The resuh 
from the study of strong gravitational lenses by Auger 
et al. (2010) is shown as a green cross in Fig. [S] Our re- 
sults are also in good agreement with those from Button 
et al. (2012) for the densest early-type galaxies in the SDSS 
(cyan dot). 

As a final remark, we note that the absence of adiabatic 
contraction indicates that non-dissipative mergers and/or 
feedback play an important role in the formation of early- 
type galaxies of all masses. 
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APPENDIX A: CONSTRAINED ACDM BASED 
MASS MODELS OF EARLY- TYPE GALAXIES 

In this appendix we describe the mass models we construct 
to reproduce the observed structural and dynamical scal- 
ing relations of early-type galaxies. These are the same as 
described in Dutton et al. (2011a). Our most general mass 
model consists of three spherical components: a stellar com- 
ponent with Sersic index n = 4; a stellar component with 
Sersic index n = 1; and a dark matter halo. The dark matter 
halo is an NEW (Navarro, Frenk, & White 1997) modified 
by the response of the halo to galaxy formation. This model 
has 7 parameters (4 for the stars and 3 for the dark mat- 
ter), 5 of which are determined by observations and theory. 
The distribution of stellar mass in galaxies is described by 
three relations: half-light size vs SPS mass (for the n = 1 
and n = 4 components) and de Vaucouleurs fraction /dcv vs 
SPS mass (the fits to these relations are given in Table 1X2)) . 
The relation between SPS mass and dark halo mass is taken 
from the compilation of observations by Dutton et al. (2010). 
The relation between dark halo concentration and dark halo 
mass is taken from ACDM N-body simulations in a WMAP5 
cosmology (Maccio et al. 2008). 

The two unknowns are the normalization of the stel- 
lar masses, which we term Aimf, and the response of the 
dark matter halo to galaxy formation. Following Dutton 
et al. (2007; 2011) we consider a range of halo responses 
parametrized by the parameter v. Where standard adiabatic 
contraction (Blumenthal et al. 1986) corresponds to v = 1, 
the contraction model of Gnedin et al. (2004) corresponds 
to ~ 0.8, the contraction model of Abadi et al. (2010) 
corresponds to ~ 0.4, an unmodified halo corresponds to 
u = Q, while expansion corresponds to < 0. As a limiting 
case of maximum halo expansion we also consider models in 
which mass follows light (i.e., no dark matter). 

The combinations of allowed Aimf and v are con- 
strained by the observed velocity dispersion vs SPS mass 
relation (Fig. [T] Table In Dutton et al. (2011a) we 
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used empirical constraints (from strong lensing and dynam- 
ics) for the average conversion between circular velocity at 
the half-light radius Vcirc(-Rc) and the velocity dispersion 
within the half-light radius, tJo. In this paper we compute 
aperture velocity dispersions for our models by solving the 
spherical Jeans equations, as describe below. 



Al Predicting SDSS aperture velocity 
dispersions. 

As described in the previous section, we have spherical mod- 
els for the distribution of total and stellar mass. Given an 
assumption of the anisotropy profile we can then solve the 
spherical Jeans equations to get the radial velocity disper- 
sion profile. We can then compute the projected velocity dis- 
persions inside the SDSS aperture. The relevant equations 
are given below. 

We consider spherical galaxy models with 3D stellar 
density distribution p*{r), and projected stellar density dis- 
tribution E*(_R). The radial component ar{r) of the velocity 
dispersion tensor is found by solving the Jeans equations: 

d(p^ 2/3 2 GMjr) 

where /3 is the velocity anisotropy, and M{r) is the spher- 
ically enclosed mass within radius r. The solution of the 
Jeans equation is given by 



2 1 



(A2) 



where I — exp J{2l3/r)dr is the integrating factor. For con- 
stant P the integrating factor is I = r'^^ If we start with 
a 2D surface density profile (e.g., a Sersic profile), the 3D 
surface density profile is given by 



dE,(i?) dR 



dR ^ 



(A3) 



Alternatively, if we start with a 3D density profile (e.g., a 
Hernquist profile), the projected density profile is 



E,(i?) 



p.(r) 



: dr. 



>R v -R^ 
The projected velocity dispersion is given by 



E.cTp = 2 



1 



4' 



p*CTr 



i?2 



: dr. 



(A4) 



(A5) 



The SDSS spectra are measured within a 3 arcsec diameter 
aperture, so the aperture projected velocity dispersion is 
given by 



^{E. a^) 2tv RdR 
l^''''{'S,)2nRdR 



(A6) 



Where the angled brackets indicate variables that have been 
convolved with the seeing. We assuming a Gaussian with 
FWHM=1.4 arcsec, which is the median seeing for SDSS 
(Simard et al. 2011). 

Fig. lAll shows an example mass model consisting of an 
NFW dark matter halo and a Hernquist bulge of stars. The 
upper panel shows circular velocity (black lines) and aper- 
ture velocity dispersion (red lines) profiles. The total circu- 
lar velocity is given by the solid black line, while the con- 
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Figure Al. Example circular velocity and projected velocity- 
dispersion profiles for a model galaxy made of a Hernquist sphere 
of stars and a cosmologically motivated NFW dark matter halo. 



tribution from the stars is given by the black dashed line. 
The dark matter halo makes up the difference between these 
two lines. The aperture velocity dispersions are given for 
different values of the anisotropy parameter, /?, chosen to 
bracket the values observed in early-type galaxies (Gerhard 
et al. 2001). 

The vertical lines mark useful apertures: -Re is the pro- 
jected half-light radius, R1/2 is the 3D half-light radius. The 
two most common apertures for measuring or "correcting" 
aperture velocity dispersions are Rc and Rc/8. 

As shown by (e.g.. Wolf et al. 2010), if the radial veloc- 
ity dispersion profile is constant (i.e., d\na^/dlnr = 0), 
then the mass within the 3D half-light radius Mi/ 2 ~ 
3(jfo3_Ri/2/G, where cr'^los is the integrated line-of-sight ve- 
locity dispersion of the system. This can be more compactly 
expressed in terms of circular velocity: V1/2 = \/3o"k)s- How- 
ever, in general galaxies do not have constant velocity disper- 
sion profiles, and for galaxies in SDSS the observed aperture 
typically contains less half the light (see below). Thus we ex- 
pect some variation from the formula from Wolf et al. (2010) . 

The lower panel in Fig. lAll shows the relation between 
the circular velocity at the 3D half-light radius (V1/2), and 
the projected velocity dispersion ((Tap(-R)) measured within 
a radius R. For large apertures R ~ 2Rc the circular veloc- 
ity is indeed almost independent of anisotropy. However, for 
the typical apertures for SDSS galaxies (0.5 ^ Rc ^ 1-0), 
there is a non-negligible dependence of circular velocity on 
anisotropy. Furthermore, when using velocity dispersions 
within on eighth of an effective radius, there is a large de- 
pendence of circular velocity on anisotropy Vi/2/o-ap varies 
from 1.3 to 1.8. Thus by correcting the aperture to smaller 
radii than observed information is lost, while correcting the 
aperture to larger radii than observed a specific mass profile 
is assumed. This motivates us to model the aperture velocity 
dispersions directly. 
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Table Al. Observed aperture velocity dispersion (o"ap) - stellar mass (Mspg) and circularized effective radius {Rc) - stellar mass relations 
from Fig. [l] and the correlation between tlie residuals of these relations from Fig. [3l 



logio(Msps/A^0) log 10 ((Tap/ km s ^) logm (flo/kpc) 

min max median median 16th % ile 84th % ile median 16th % ile 84th % ile 



Alogio(o-ap)/Alogio(_Re) 



9.4 
9.6 
9.8 
10.0 
10.2 
10.4 
10.6 
10.8 
11.0 
11.2 
11.4 
11.6 



9.6 
9.8 
10.0 
10.2 
10.4 
10.6 
10.8 
11.0 
11.2 
11.4 
11.6 
11.8 



9.515 
9.712 
9.912 
10.112 
10.313 
10.508 
10.707 
10.903 
11.098 
11.293 
11.486 
11.678 



1.748±0.009 
1.822±0.005 
1.896±0.004 
1.977±0.002 
2.056±0.001 
2.119±0.001 
2.179±0.001 
2.238±0.001 
2.294±0.001 
2.347±0.001 
2.392±0.001 
2.431±0.001 



1.602 
1.710 
1.786 
1.862 
1.955 
2.031 
2.093 
2.155 
2.218 
2.280 
2.333 
2.375 



1.856 
1.942 
2.028 
2.091 
2.153 
2.208 
2.261 
2.310 
2.359 
2.405 
2.447 
2.482 



0.078±0.012 
0.083±0.009 
0.139±0.006 
0.199±0.004 
0.260±0.002 
0.348±0.002 
0.453±0.001 
0.564±0.001 
0.676±0.001 
0.791±0.001 
0.921±0.001 
1.064±0.002 



-0.133 
-0.106 
-0.051 
0.012 
0.106 
0.206 
0.324 
0.450 
0.568 
0.688 
0.814 
0.951 



0.236 
0.256 
0.313 
0.361 
0.425 
0.498 
0.592 
0.691 
0.796 
0.908 
1.042 
1.181 



-0.368±0, 
-0.385±0, 
-0.443±0, 
-0.406±0, 
-0.377±0, 
-0.357±0, 
-0.337±0, 
-0.304±0, 
-0.261±0, 
-0.197±0, 
-0.152±0, 
-0.148±0, 



018(lHf4) 

oi2(l°:°^i) 

009(1^1?) 

oo4(i°:»i) 
oo4(i°:»^) 

oo3(l°:»?) 
oo3(ig:»i) 



Table A2. Parameters of the fitting formula (Eqs. [4]& Ell for various scaling relations used in this paper. 
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XQ 


yo 
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logio(a-ap/kms-l) 


logio(MsPs/M0) 


0.42 


0.20 


10.7 


2.18 


1.0 


logio(iJo/kpc) 


logio(A^SPs/M0) 


0.00 


0.70 


10.2 


0.22 


1.0 


logioC-f^ap/kpc) 


logio(Msps /Mq) 


0.45 


0.24 


11.1 


0.54 


1.2 


'^log flap 


logio(Msps/M0) 


0.046 


-0.037 


10.1 


0.098 


2.0 


/doV 


logio(MsPs/M0) 


0.17 


-0.22 


10.9 


0.65 


2.0 


logio(-Rdcv/kpc) 


logio(Msps /Mq) 


-0.01 


0.55 


10.0 


-0.08 


2.0 


logio(-Rcxp/kpc) 


logio(MsPs /Mq) 


0.29 


0.61 


10.1 


0.39 


1.0 


e log (Tap 


logio(o-ap/kms-i) 


-0.8 


0.0 


1.7 


0.093 


3.8 


e log _Rc 


logio{_Rc/kpc) 








0.035 




e log Ms PS 


loglo(MsPS /Mq) 








0.092 
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yi 
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7 


""logioCo-ap) 


logio(A^SPs/M0) 


0.12 


0.045 


10.7 




1.0 


°'loglo(«c) 


logio(A^SPs/M0) 


0.18 


0.10 


10.55 




1.8 



A2 Physical aperture sizes 

Since we are modeling the aperture velocity dispersion, we 
need to know the physical size of the aperture. This will de- 
pend on the redshift of the galaxy. Fig. IA2I shows the aper- 
ture radius in kpc as a function of stellar mass. In the upper 
panel the solid points show the median relation, while the 
solid line is a fit to these data points (see Table IA2I for the 
parameters of the fit). The red dashed line shows the me- 
dian effective radius (from Fig. which shows that for all 
but the lowest and highest mass galaxies the average SDSS 
aperture is close to the half-light radius. For the most mas- 
sive galaxies the median aperture is roughly 0.5-Rc- There is 
scatter in the aperture radius at fixed mass, but it turns out 
not to induce any significant scatter to the model velocity 
dispersions. 

A3 Sampling of Galaxies 

The observed sampling of the stellar masses of galaxies is 
not uniform. The distribution is very roughly log-normal 
with a peak at Afsps = 10^^ Mq (error bars in Fig. IA3|) . 



Measurement uncertainties in stellar masses will preferen- 
tially scatter galaxies away from the peak. We fit the distri- 
bution of stellar masses with a set of Gaussians (red lines 
in Fig. IA3|I , and since the measurement errors are roughly 
constant, we can thus deconvolve the observed distribution 
of stellar masses analytically (green lines in Fig. IA3[) . The 
differences between the model and deconvolved model are 
shown in the lower panel of Fig. IA3I The most significant 
differences occur for high masses (Msps > 10"-'^ Mq), where 
the deconvolved model has up to a factor of 1.6 fewer galax- 
ies. The net effect of the non-uniform stellar mass distribu- 
tion coupled to errors on stellar masses is to make the ob- 
served relations (slightly) shallower than the true relations. 



A4 Adding scatter to the models 

We add both intrinsic and observational sources of scatter 
to our models. There are three primary sources of intrinsic 
scatter: 1) Scatter in galaxy size at fixed stellar mass; 2) 
Scatter in the "pristine" dark halo concentration at fixed 
halo mass; and 3) Scatter in dark halo mass at fixed stellar 
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Figure A2. Aperture size - stellar mass relation. Upper panel: 
Solid points show median of aperture sizes in kpc, solid line shows 
a fit to these points. The dotted lines show the 16th and 84th 
percentiles. For comparison the red dashed line shows the median 
half-light radius. This shows that for all but the highest and lowest 
masses, the SDSS fiber aperture, on average, corresponds to the 
half-light radius. Lower panel: Scatter in the aperture sizes. 
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Figure A3. Observed distribution of galaxy stellar masses (as- 
suming a Chabrier IMF). Upper panel: black points with error 
bars show the data, red lines show a triple Gaussian fit, green 
lines show the intrinsic distribution after correcting for measure- 
ment errors of 0.1 dex in stellar mass. Middle panel: Same as 
upper panel but for a logarithmic scale. Lower panel: Difference 
between logjQ N of the model and deconvolved model. At high 
stellar masses the differences can be up to a factor of ^ 1.6. 



Figure A4. Measurement errors on velocity dispersion (fTap), 
half-light size (-Re), and stellar mass (Mgpg). Typical errors are 
0.03 dex on velocity dispersion, 0.03 dex on size, and 0.1 dex 
on stellar mass. The average errors are fitted with Eq. |4l with 
parameters given in Table |A2] 

mass. The scatter in size at fixed stellar mass is taken from 
the observed RM relation (Fig. [T|. The scatter in the pris- 
tine halo concentration for relaxed haloes in cosmological 
simulations is ~ 0.11 dex (Jing 2000; Wechsler et al. 2002; 
Maccio et al. 2008) . The scatter in stellar mass at fixed halo 
mass is about 0.15 dex (e.g., More et al. 2011). However 
due to shallow slope of the stellar mass vs halo mass rela- 
tion at high halo masses, the scatter in halo mass at fixed 
stellar mass increases with stellar mass. The results from 
More et al. (2011) for red galaxies can be approximated by 
flog A/^ir = 0.15 -f 0.2(log Msps — 10.4), with a minimum of 
0.15. ™ 

Additionally there is scatter in the aperture radius at 
fixed stellar mass from the redshift sampling of galaxies, and 
is given in Fig. IA2I It turns out that this adds negligible 
scatter to the VM relation, and so we do not include it in 
our fiducial models. 

There are measurement errors in stellar masses, sizes 
and velocity dispersions. These measurement errors are 
roughly constant (Fig. IA4( ) with aiogMgpg ~ 0.1, oiosR^ — 
0.035 and fiogo-j^p ~ 0.04. For our purposes the measure- 
ment errors in stellar masses and sizes are more important 
than the measurement errors in velocity dispersions. This is 
because we bin galaxies in both stellar mass and size, but 
not velocity dispersion. Errors in stellar masses can change 
the slope of the correlation between VM and RM residuals, 
and thus are especially important to account for. 
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